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ABSTRACT 

We report the detection of three new extrasolar planets orbiting the solar type stars HD 85390, HD 90156 and HD 103197 with the 
HARPS spectrograph mounted on the ESO 3.6-m telescope at La Silla observatory. HD 85390 has a planetary companion with a 
projected intermediate mass (42.0 M B ) on a 788-day orbit (a=1.52 AU) with an eccentricity of 0.41, for which there is no analogue in 
the solar system. A drift in the data indicates the presence of another companion on a long period orbit, which is however not covered 
by our measurements. HD 90156 is orbited by a warm Neptune analogue with a minimum mass of 17.98 M e (1.05 Mg), a period 
of 49.8 days (a=0.25 AU) and an eccentricity of 0.31. HD 103197 has an intermediate mass planet on a circular orbit (P = 47.8 d, 
M sin i = 3 1 .2 M e ). We discuss the formation of planets of intermediate mass (~ 30 - 100 M B ) which should be rare inside a few AU 
according to core accretion formation models. 

Key words, stars: planetary systems - techniques: radial velocities - stars: individual: HD 85390 - stars: individual: HD 90156 - 
stars: individual: HD 103197 - stars: planetary systems: formation 



1. Introduction 

The HARPS (High Accuracy Radial Velocity Planet Searcher) 
program has been on-going since 2003 at the ESO 3.6 meter 
telescope located at La Silla Observatory in Chile (Pepe 2002). 
Thanks to its high accuracy of < 1 m/s and its long time sta- 
bility, HARPS has in particular detected many low mass planets, 
starting with fi Arae b at the beginning of the program (Santos 
et al. 120041 1 to the recent two Earth-mass planet around GJ 581 
(Mayor et al. 2009). As such low mass planets seem to be very 
abundant (Lovis et al. [2009), HARPS helps significantly to in- 
crease the number of known extrasolar planets. The latter is rais- 
ing rapidly in time f, roughly as f 2 3 since the discovery of 51 
Peg b by Mayor & Queloz ( 1995 ) fifteen years ago. 

These detections of extrasolar planets have enormously stim- 
ulated the research on planet formation. They have shown the 
importance of mechanisms which were previously underesti- 
mated in their importance from the shape of our own solar 



* Based on observations made with the HARPS instrument on the 
ESO 3.6 m telescope at La Silla Observatory under the GTO programme 
ID 072.C-0488. 

** Individual radial velocities are available electronically at CDS via 
anonymousftptocdsarc.u-strasbg.fr (138.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/ 



system alone, like for example orbital migration which was 
known on purely theoretical bases for along time (Goldreich & 
Tremaine, 1980), overthrowing the idea that giant planets can 
only be found beyond a few astronomical units (Boss 1995). 

The quickly growing number of extrasolar planets has al- 
lowed to derive many important distributions of physical proper- 
ties and orbital elements. We know now distributions of masses, 
semi-major axes or eccentricities of giant planets within a few 
AU around solar like stars (e.g. Udry & Santos 2007 ), and see 
correlations between stellar and planetary parameters. The most 
important example is the "metallicity effect" i.e. the increase of 
giant planet frequency with stellar [Fe/H], see e.g. Santos et al. 
(EOTJfl or Fischer & Valenti ([20051) . 

By varying initial conditions of planetary formation mod- 
els in a Monte Carlo way (Ida & Lin 2004 ), one can try to de- 
rive these statistical distributions also from theoretical forma- 
tion models and compare them with the observed ones. This 
young method of planetary population synthesis (Mordasini et 
al. [2009a, 2009b) allows to use the full wealth of observational 
data to constrain theoretical models and to improve our under- 
standing of planet formation, as demonstrated here for the spe- 
cific case of intermediate mass planets. 

In this paper we report the discovery of two types of planets, 
namely of one Neptunian planet and of two intermediate mass 
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Table 1. Observed and inferred stellar parameters for the three 
host stars. Photometric and astrometric quantities were taken 
from the Hipparcos catalogue (ESA 1997) while stellar physical 
quantities come from Sousa et al. (2008), except for HD 103197. 
For this star, the values are from the present paper, following the 
analysis in Sousa et al. (2008 ). 



Parameter 


Unit 


HD 85390 


HD 90156 


HD 103197 


Spectral type 




K1V 


G5V 


KlVp 


V 


[mag] 


8.54 


6.92 


9.40 


B-V 


[mag] 


0.86 


0.66 


0.86 


71 


[mas] 


29.45±0.84 


45.26±0.75 


20.27±1.47 


My 


[mag] 


5.89 


5.20 


5.77 




[K] 


5186±54 


5599±12 


5303±58 


l°g£ 


[cgs] 


4.41±0.09 


4.48+0.02 


4.40±0.11 


[Fe/H] 


[dex] 


-0.07±0.03 


-0.24±0.01 


0.21±0.04 


L 


[L ] 


0.43±0.03 


0.72±0.01 


0.47±0.11 


M„ 


[M e ] 


0.76 


0.84 


0.90 


v sin 


[km/s] 


1 


<1 


2 


log R HK 




-4.97 


-4.95 


-5.07 


P mt 


[days] 


44±5 


26±3 


51+5 


Age 


[Gyr] 


~ 7.2 


~ 4.4 


~ 9.1 



planets (about 31 and 42 M m ) without analogue in the solar sys- 
tem. We discuss the formation of the latter objects in the context 
of the core accretion formation model at the end of the paper. 
Such intermediate mass planets are interesting, as they provide 
observational constraints on the rate at which gas can be accreted 
by planets starting runaway gas accretion. 

The paper is organized as follows: In section|2]we discuss the 
host star properties. Radial velocity measurements and orbital 
solutions are presented in section [3] In section [4] we summarize 
our discoveries and in section[5]we discuss them. 



2. Host star characteristics 

Basic photometric and astrometric properties for the three 
host stars HD 85390 (HIP 48235), HD 90156 (HIP 50921) and 
HD 103197 (HIP 57931) are taken from the Hipparcos catalogue 
(ESA 1997). Accurate spectroscopic stellar parameters and de- 
rived quantities for HD 85390 and HD 90156 are taken from 
the analysis of Sousa et al. (2008). Sousa et al. (2008) used 
high quality, high signal to noise HARPS spectra to obtain ho- 
mogenous estimates of the spectroscopic stellar parameters of all 
451 targets in the HARPS Guaranteed Time Observations (GTO) 
"high precision" sample, to which the two mentioned stars be- 
long. HD 103197 is in contrast part of the lower RV precision 
"CORALIE extension" sample, see Naef et al. (2007). For this 
star, we stacked several individual spectra together to obtain a 
high signal to noise spectrum and repeated the analysis described 
in Sousa et al. (20Q8]>. 

Individual HARPS spectra were used to derive the radial ve- 
locity RV, the Bisector Inverse Slope of the Cross-Correlation 
Function, the vsin/ of the star as well as the chromospheric 
activity index \ogR HK , using a similar recipe as Santos et al. 
(2000). From the activity indicator we also derive an estimate 
of the stellar rotation period P rot using the empirical Noyes et 
al. (119841) activity-rotation correlation. Approximate ages are 
obtained using the improved activity-rotation-age calibration of 
Mamajek & Hillenbrand (2008). These ages typically have er- 
rors of about 0.2 dex. 



The basic properties of the four stars are summarized in 
Table[T] The stellar mass estimates from Sousa et al. (2008|) typ- 
ically have an error of 0.1 M Q . 



2.1. HD 85390 (HIP 48235) 

HD 85390 is an Kl dwarf in the southern hemisphere. The 
Hipparcos catalog (ESA 1997) lists a visual magnitude V = 
8.54, a color index B - V — 0.86, and an astrometric parallax 
n = 29.45 ± 0.84 mas, setting the star at a distance of 34 ± 1 
pc from the Sun, a typical distance for planet host stars in RV 
search programs. Its absolute magnitude is then estimated to be 
My = 5.89 mag. The spectral analysis of Sousa et al. ([2008) 
yields the following physical parameters: An effective tempera- 
ture r e ff = 5186 ± 54 K, a surface gravity logg = 4.41 ± 0.09, 
and a metallicity [Fe/H]= -0.07 + 0.03 dex. The star has thus a 
metallicity similar to the Sun. The derived mass is 0.76 M , and 
the luminosity is estimated from the absolute magnitude, the ef- 
fective temperature and the corresponding bolometric correction 
to be 0.43 Lq, thus less than half the solar luminosity. A low pro- 
jected rotational velocity vsin / of about 1 km/s is derived from 
a calibration of the width of the cross-correlation function as de- 
scribed in Santos et al. ((2002 ). With an activity indicator log R HK 
of -4.97, the star is among the non-active stars in our sample, and 
no large radial velocity jitter is expected. The detailled way how 



we deal with the stellar jitter is described in section 3.1 

Several additional measurements of the properties of 
HD 85390 can be found in the literature that agree fairly well 
with the mentioned values: Gray et al. (2006 ) list a spectral type 
K1.5V, r eff = 5069 K, logg = 4.48, [M/H]= -0.11 and \ogR HK 
of -5. 06. He nry et al. (119961) measure a log R HK of -4.93. Minniti 
et al. ( 120091 ) classify the star as a Kl radial velocity stable star at 
a precision of about 5 m/s. 



2.2. HD 901 56 (HIP 50921) 

According to the Hipparcos catalogue (ESA 1997) is HD 90156 
a quite close G5V star with an astrometric parallax of n = 
45.26 ± 0.75 mas with basic photometric properties of V = 6.92 
and B - V=0.66. Sousa et al. (120081 derived from the HARPS 
spectra a T eff = 5599 + 12 K, logg = 4.48 + 0.02, and a subsolar 
metallicity of [Fe/H] = -0.24 + 0.01 dex which is well known to 
play an important role in planet formation. Other derived quan- 
tities are M* = 0.84 M Q and L = 0.72 ± 0.01 L e . The star has a 
low projected rotation velocity of less than 1 km/s which is thus 
difficult to determine accurately, see Melo et al. (2001 ) and a low 



activity index of log R , 



-4.95. 



The star has been studied by several other groups: It was 
in particular spectroscopically analyzed by Valenti & Fischer 
(2005 ). By fitting synthetic spectra they find an effective tem- 
perature of 5626 + 44 K and logg = 4.63 + 0.06. This sur- 
face gravity is significantly larger than our value or their value 
calculated from isochrone interpolation (logg = 4.46), corre- 
sponding to a systematic discrepancy noted by the authors them- 
selves. The mass derived from their spectroscopic gravities is 
with 1.25 + 0.17 M correspondingly considerably larger than 
our value. On the other hand is the mass they derive from in- 
terpolating isochrones with 0.90 ± 0.04 M similar to our result. 
Valenti & Fischer (2005) also measure a subsolar metal content 
of [Fe/H]=-0.21 ± 0.03 and estimate the age of the star to be 7.8 
Gyr. 

The metal poor nature of the star is further confirmed by 
Neves et al. ([2009) who note that the star is a member of the 
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thin disk and by Gray et al. (2006) who find a [M/H] of -0.26. 
The latter authors also give a r e ff = 5578 K, logg = 4.52 and 
\ogR HK = -4.96 which is comparable to our values. HD 90156 
is listed in the catalogue of suspected variable stars of Kukarkin 
et al. ( 119811 1 with an amplitude in V of 0. 1 mag. The much more 
accurate Hipparcos data does however not point to any signifi- 
cant variability. 

2.3. HD 1031 97 (HIP 57931) 

HD 103197 is classified in the Hipparcos catalogue (ESA 1997) 
as KlVp, with V=9.40, B-V=0.86 and n = 20.27 + 1.47 mas. 
Repeating the same spectroscopic analysis as described in Sousa 
et al. d20Q81 > with stacked HARPS spectra leads to the following 
values: r eff = 5303 ± 58 K, logg= 4.40 + 0.1 1 and a supersolar 
metallicity of [Fe/H]= 0.21 + 0.04 dex. The derived mass and 
luminosity is 0.90 M Q and 0.47 L . From the HARPS spectra we 
also infer the following quantities: a projected rotational velocity 
of roughly v sin i « 2 km/s, a low activity indicator log R HK = 
-5.07 and an estimated P m of about 51+5 days. This star is 
therefore very inactive. 

For comparison, Jenkins et al. (2008) find also a low 
log R HK = -5.05, but their [Fe/H]=-0.12 ± 0.04 is clearly lower 
than our value. 



3. Radial velocity measurements and orbital 
solutions 

3.1. Jitter and error estimates 

The quantified error of an individual radial velocity measure- 
ment includes photon noise, calibration and instrumental drift 
uncertainty. Additional, not quantified errors which are not in- 
cluded in this number will originate from intrinsic radial veloc- 
ity variability of the star (jitter) due to granulation (Kjeldsen et 
al. [2005 ), magnetic activity and pulsations (see Pepe & Lovis 
2008 for an overview). This jitter can in principle be a serious 
threat to high precision RV measurement at a level of 1 m/s 
(Wright 2005 ). On the other hand, there are now severals stars in 
the HARPS high precision program like for example HD 69830 
(Lovis et al. 12^061 ) or HD 40307 (Mayor et al. 120091 for which 
residuals around the fit clearly below 1 m/s have been found, like 
0.64 m/s for HD 69830 (Lovis et al. 12006] ). 

Neglecting the RV jitter in the total error estimate can have 
negative effects in the derivation of the uncertainties of the or- 
bital parameters of the fits. Neglecting it has also the effect that 
measurements with a small quantified error (particularly high 
S/N, due to e.g. very good seeing to which HARPS is sensitive) 
get an inappropriately high weight in the fitting procedure, caus- 
ing also misleading results. 

Several authors have tried to derived empirical relations link- 
ing a star's radial velocity jitter to other quantities like B - V or 
R HK (e.g. Saar et al. 119981 Wright et al. 120051 ). Concerning the 
latter work, all three stars presented here would fall in the class 
with the lowest derived median jitter of 3.5 m/s. The jitter of the 
stars presented here is however obviously significantly smaller, 
as can be seen by the fact that planets which induce a radial ve- 
locity semi-amplitude of that order are clearly detected in the 
data. This indicates that those older jitter estimates were clearly 
affected by instrumental and data reduction uncertainties. 

We find from our HARPS measurements that for the typical 
stars in the high precision program with insignificant evolution, 
a very low activity level (log R HK < -4.95 like for the three stars 



in this work), a small intervall in B - V (usually 0.5-1.0, 0.66- 
0.86 in this paper), and a small v sin i (< 2 km/s), it is difficult to 
derive empirical relations like the ones mentioned before which 
cover much larger domains of these quantities. For example, we 
do not see an obvious correlation between the minimal RV jitter 
and log R HK or B - V among the stars with the mentioned prop- 
erties in our sample. In a recent study, Isaacson & Fischer (2010) 
found at least qualitatively very similar conclusions, with in par- 
ticular a jitter level independent of activity indicators and below 
their instrumental and reduction precision (~ 1.6 m/s) for mid 
to late K dwarfs. For somewhat bluer stars, they derived slightly 
higher astrophysical jitter levels (~ 1 .4 m/s), but still only a weak 
correlation with chromospheric activity. 

It seems that in the HARPS high precision data, the lower 
boundary of the minimal total rms is of the order of 50 cm/s 
for stars with parameters mentioned above, and our current re- 
duction pipeline. As gaining insight in the fundamental limiting 
factors for the radial velocity technique is an important task, we 
will present an in-depth analysis of the statistical properties of 
the sample regarding RV jitter in a dedicated forthcoming paper 
(Lovis et al. in prep.). 

For the work presented here, we restrict ourselves to a prag- 
matic approach and compare the three stars with similar ones in 
the sample for which the jitter can be estimated in form of the 
rms of the residuals to the fitted orbits, cr(O-C) (this number can 
however also contain additional contributions e.g. from unde- 
tected low mass planets). Both HD 85390 (Kl V) and HD 103197 
(K1V) are similar to HD 69830 (K0V) in terms of spectral type, 
effective temperature and logR HK , even though that HD 85390 
is somewhat cooler than this star. As mentioned, for HD 69830 a 
<j{0-C) of 0.64 m/s was found in Lovis et al. (2006 ) for the later 
measurements. Note that since then, our calibration algorithm 
has been improved even more (Pepe & Lovis 120081 . HD 90156 
(G5V) is similar to HD 47 186 (G5V) for which a total cr{0 - C) 
of 0.91 m/s was found in Bouchy et al. (2009). From this we 
conclude that the typical levels of the unquantified errors are of 
order of 0.5 to 0.8 m/s, where the latter value was already used 
in Pepe et al. d2007l . 

Following Saar et al. (11998b or Pepe et al. (120071 ) we then as- 
sume that we can quadratically add this estimated radial velocity 
jitter to the quantified error of each radial velocity measurement, 
using a constant value of either 0.5 or 0.8 m/s. This might not be 
completely justified as depending on the physical nature of the 
variability source like magnetic cycles of the star (Santos et al. 
2010), the jitter might be rather of a red noise type, and vary- 
ing in time. Using such an additive procedure, we can however 
certainly mitigate the adverse effects of artificially high dynam- 
ics in the quantified (photonic) errors. We then search for orbital 
solution in the data after adding the two mentioned values, and 
compare the impact on the best fit parameters and their errors. 

3.2. HD 85390: An intermediate mass planet and a long-term 
drift 

HD 85390 was monitored with HARPS for a time span of 6.5 
years (2386 days), during which 58 measurements were ob- 
tained. Except for the first six measurements we always fol- 
lowed the observational strategy to use a fixed exposure time 
of 15 minutes in order to minimize the effects of stellar acoustic 
modes (Pepe & Lovis |2008). We however nevertheless include 
the early, less accurate measurements to have a better constraint 
on the long time drift present in the data, cf. next. The spectra 
which were all obtained in ThAr simultaneous mode have typical 
S/N at A — 550 nm of 1 16 and a mean quantified error of the ra- 
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dial velocity (including photon noise, calibration and instrumen- 
tal drift uncertainty) of 0.62 m/s. Given the low chromospheric 
activity of the star, they should however induce a clearly smaller 
raw rms than the one found in the 58 measurements which is 3.3 
m/s (peak-to-peak 14.4 m/s), indicating instead the presence of 
a companion. To account for the unquantified errors (see section 
|3.1| >, we add quadratically to all our measurements an additional 
error of 0.5 resp. 0.8 m/s. Using either 0.5 m/s or 0.8 m/s has 
however only very minor influences on the orbital solutions: the 
derived quantities like the period P, the mean longitude A, the 
eccentricity e or the velocity semi-amplitude K agree to better 
than 2 % for the two cases, much less than the errors bars on 
these quantities. This is in good agreement with Wright (2005) 
who states that radial velocity jitter is more important in the de- 
tection of extrasolar planet than in their characterization, and that 
the effects of including jitter in the noise estimation for the de- 
rived orbital parameters of the best fit is typically small, at least 
for a case like here where the induced radial velocity signal is 
larger than the estimated jitter amplitude. 

Radial velocity measurements as a function of Julian date 
(with 0.8 m/s of additional error) are shown as red dots in the 
top panel of Fig. T] A long-term, linear drift overlaid by a higher 
frequency modulation with a period of about 800 d can be seen 
directly in the data. Of the latter, three phases are covered. The 
bottom panel of Fig. T] shows the generalized Lomb Scargle pe- 
riodogram (as described in Zechmeister & Kiirster 2009) of the 
data. The peak at such a period is clear and has a very low FAP 
of less than 0.01%, as found by scrambling the data 10000 times 
which confirms the visual impression of a very clearly detected 
orbit. The FAPs are obtained by performing random permuta- 
tions of the velocities, calculating the periodogram, recording 
the peak power for each trial, and finally comparing the power 
of the real signal to the peak power distribution of the permuted 
datasets. 

Therefore we fit the measurements with a linear drift and a 
Keplerian orbit, using a genetic algorithm (see below for a dis- 
cussion of the nature of the drift). The resulting best fit and the 
residuals around it are also plotted in the top panel of Fig. [T] 
The phase-folded radial velocity curve of the planet candidate 
after removing the drift is shown in the middle panel of this fig- 
ure. The reduced^ per degree of freedom is 1.52 (2.31 for an 
additional error of 0.5 m/s) and the weighted r.m.s. of the resid- 
uals around the solution cr(0 - C) is 1.15 m/s (1.14 m/s for an 
additional error of 0.5 m/s), which is not much more than the 
mean total assumed error of the individual measurements. In the 
residuals, no additional significant periodic signals is found. 

The somewhat higher cr(O-C) is also at least partially due to 
the mentioned six early measurements which can be identified in 
Fig. [T]by their larger error bars. Excluding these six early radial 
velocities results in a very similar fit with a reduced^ = 1.10 
(additional error of 80 cm/s) and reduces cr(0 - C) to just 1.00 
m/s. 

The long term signal is found to have slope of 1 .15+0.02 m/s 
yr _I . The derived parameters for the Keplerian orbit imply a pe- 
riod P = 788+25 d (corresponding to a semi-major axis a — 1.52 
AU), a minimum mass of Msin; = 42.0 + 3.6 M e (2.44 Mg), 
and eccentricity e = 0.41 + 0.12. The error intervals were com- 
puted for a 68% confidence level with 5000 Montecarlo itera- 
tions, where the actual measurements are replaced with a value 
drawn from a normal distribution around the nominal value with 
a standard deviation which is the same as the error bar on the 
point. These altered data sets are then again fitted to obtain the 
confidence intervals from the distribution of orbital elements one 
finds. 
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Fig. 1. Radial velocity measurements of HD 85390 as a func- 
tion of Julian Date (top). The six earliest measurements have 
larger error bars. The best solution with a linear drift plus one 
Keplerian orbit is also displayed as well as the residuals to this 
model. The middle panel shows the radial velocity as a func- 
tion of the orbital phase for HD 85390b, with the effect of the 
drift removed. The bottom panel shows the generalized Lomb 
Scargle periodogram of the observations. The peak at a period 
of about 800 days is very prominent, and has a false alarm prob- 
ability (FAP) of much less than 0.1 %. The three horizontal lines 
indicate FAPs of 0.1, 1 and 10% from top to bottom. 



The orbital elements are listed in Table [2] We checked that 
there is no correlation between bisector shape of the cross- 
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Table 2. Orbital and physical parameters of the 1 -planet plus 
linear drift solution of the HARPS data for the objects around 
HD 85390. A is the mean longitude at the barycenter of the ob- 
servations at BJD 54269.3584. 



HD85390 



HARPS+CORALIE 



Parameters 




HD 85390 b 


P 


[days] 


788 ± 25 


A 


[deg] 


177 ± 25 


e 




0.41 ±0.12 


CO 


[deg] 


-94 ± 23 


K 


[m/s] 


3.82 ±0.33 


V 


[km/s] 


33.0853 ± 0.0004 


slope 


[m/s yr 1 ] 


1.15 ±0.02 


f(m) 


[i<r 9 M ] 


0.0035 ± 0.0009 


M sin i 


[M ffi ] 


42.0 ±3.6 


a 


[AU] 


1.52 ±0.04 


N 




58 


span 


[years] 


6.5 


cr(0 - C) 


[m/s] 


1.15 






1.52 



correlation function and the radial velocity or the residuals of 
the Keplerian fit, providing strong support to the planetary inter- 
pretation of the RV signal. 

HD 85390 is also p art of the CORALIE planet search sur- 
vey (Udry et al. 120001 ). Including the 19 measurements made 
with CORALIE expands the observed time base significantly to 
11.3 years (41 14 days), back to JD=2 451 262 which could help 
to better constrain the long-term drift seen in the HARPS data. 
The radial velocity as a function of time for the two instruments 
is shown in Fig. [2] It is found that at the level of precision of 
CORALIE (about 6.5 m/s for this target), the star is stable, i.e. 
the ratio of the rms of all the velocity measurement and the mean 
error of one velocity measurement is about unity. Repeating the 
fitting process with the combined data set leads to a solution for 
the Keplerian orbit which agrees within the error bars with the 
solution without the CORALIE data, namely P = 821 ± 20 d, 
e = 0.34 ± 0.10 and Msin/ = 42.3 + 3.5 M m , with a x 2 r of 1.59 
for additional jitter of 0.8 m/s in the HARPS measurement. The 
linear component of the drift is now 1.03 + 0.02 m/s yr -1 . From 
that we can exclude the presence of a massive outer companion 
with a period smaller than about 10 years. 

In the nominal model for the orbital solution, we have as- 
sumed a Keplerian orbit plus a linear drift. One wonders if we 
could significantly improve the fit by allowing for a quadratic 
component in the drift, as one could tentatively infer from two 
early CORALIE measurement lying 12 resp. 16 m/s above the 
linear fit. 

This could be indicating the presence of a body that might 
be interesting to direct imaging searches. To explore this sig- 
nificance, we have fitted the combined dataset of HARPS and 
CORALIE which is more relevant in this context as it constrains 
better the long time behaviour also with a Keplerian plus a 
quadratic drift. The orbital solution for the Keplerian orbit cor- 
responds then to a companion with properties that again agree 
within the error bars with the solution shown in table |2] This 
is due to the fact that the HARPS measurements alone constrain 
this orbit. The reduced x 2 r for the combined data are now 1 .47 
for an additional errors of 0.8 m/s, i.e. somewhat less than for 
the linear model with one degree of freedom less. We have then 
used an F-test as in Pourbaix & Arenou (2001 ) to get the statisti- 
cal probability that the improvement due to the quadratic model 
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Fig. 2. Radial velocity measurements of HD 85390 as a function 
of Julian Date for the combined data set of CORALIE (squares 
for measurement before its upgrade in June 2007, and diamonds 
after this moment) and HARPS (circles), together with a one 
Keplerian plus linear drift fit (solid line). The small panel shows 
the residuals around the fit. 

is significant. We find a probability of 98.7 %, corresponding 
to about 2.5 sigma. The main assumption behind the F-Test is 
however that the model is linear and that the error distribution is 
Gaussian, which is both not strictly fulfilled, therefore this num- 
ber must be regarded with caution. Thus, and because the orbital 
solution of HD 85390b itself is not significantly affected by the 
choice, we prefer to stick for the moment to the simpler linear 
model and conclude that to constrain further the nature of the 
long term drift, more high precision observations with HARPS 
are needed. 

Long term, small amplitude signals could in principle also 
be due to magnetic cycles of the host star. While we cannot ab- 
solutely exclude this possibility, we think that it is unlikely in 
the present case, due to the following reasons: First, during the 
observed time span, log R HK decreases linearly from about -4.93 
to -5.00. Thus, the change of logR HK and the long term drift 
in the radial velocity are clearly anti-correlated. This is the con- 
trary of what is expected, as in active zones, convection tends to 
be inhibited, reducing the convective blue shift, so that a shift of 
the radial velocity into the red should occur (Santos et al. l2.0101 >. 
Second, the change of logR HK , about 0.07 dex is less than the 
value of about 0. 1 dex where we start to see clear correlations of 
the activity level and the radial velocity in the HARPS database 
(Lovis et al., in preparation). Third, the amplitude of the drift is 
large (~ 8 m/s), much more than expected for the small change 
in log R HK for a K dwarf. 

3.3. HD90156: A warm Neptune analogue 

Since the mentioned change of strategy to minimize stellar noise 
(fixed exposure times of 15 minutes), we have gathered 66 ra- 
dial velocities in ThAr simultaneous mode during a time period 
of 1614 days (4.4 years). HD 90156 is a bright star (V=6.92), 
therefore the measurements have a high mean S/N of 246 at 550 
nm, with variations of the S/N between 121 and 342, depending 
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Table 3. Orbital and physical parameter of the 1 -planet 
Keplerian solution for the planet around HD 90156. The mean 
longitude A is given at the barycenter of the observations, at JDB 
54772.0784. 



Parameters 




HD 90156 b 


P 


[days] 


49.77 ± 0.07 


A 


[cleg] 


91 ± 15 


e 




0.31 ±0.10 


(O 


[deg] 


1 13 ± 18 


K 


[m/s] 


3.69 ±0.25 


V 


[km/s] 


27.0273 ± 0.0003 


Km) 


[10- 9 M o ] 


0.00022 ± 0.00005 


M sin ;' 


[M s ] 


17.98 ± 1.46 


a 


[AU] 


0.250 ± 0.004 


N 




66 


span 


[years] 


4.4 


cr(0 - C) 


[m/s] 


1.23 






2.03 



mostly on seeing conditions. The derived mean quantified er- 
ror (photon noise, calibration and drift uncertainty) is 0.41 m/s. 
The observed raw rms for the 66 measurements is 2.6 m/s, again 
more than expected for the low stellar activity level. As before, 
we add quadratically to all measurements an additional error of 
either 0.5 or 0.8 m/s and check for the consequences. 

A periodogram of the data which is displayed in the top 
panel of Fig. [3] shows a very strong peak at a period of about 
50 days with a FAP of much less than 0.1% as obtained by 
executing 10 000 data scrambling experiments. A number of 
aliases can also been seen. A one Keplerian fit assuming an addi- 
tional RV error of 0.8 m/s quickly converges on a solution with 
a Msin/ = 17.98 + 1.46 M m planet on an orbit with a period 
of 49.77 ± 0.07 days and an eccentricity of 0.31 + 0.10. This 
solution and the measurements are shown in the lower two pan- 
els of Fig. pi The reduced x 1 ,- is 2.03, and the cr(0 - C) is 1.23 
m/s (see Table [3J. The mentioned period corresponds to a semi- 
major axis of 0.25 AU, and the projected mass is very close to the 
mass of Neptune (1.05 Mg), making this planet a warm Neptune 
analogue. We checked that there is no correlation between the 
bisector shape of the cross-correlation and the radial velocity or 
the residuals around the Keplerian fit. Assuming an additional 
error of 0.5 instead of 0.8 m/s leads to virtually identical fit pa- 
rameters with relative differences of less than 3 %. The reduced 

increases to 3.91, indicating that we now underestimate the 
unquantified errors in the measurements. We note the subsolar 
metallicity of the parent star, so that HD 90156 b is another ex- 
ample of a Neptunian planet not following the "metallicity ef- 
fect" seen for giant planets, as discussed for example in Sousa 
et al. J2008I) . It is also an example of a warm low mass planet 
which does not show the same strong pile-up behaviour at about 
0.04 AU as higher mass Hot Jupiters do (Udry & Santos 120071 

In a periodogram of the residuals around the one planet fit, 
there is a signal at a period of about 14 days. With a FAP of 16 
% we can however not yet conclude anything about its origin. 
We however hope to soon be able to either confirm or reject this 
signal by additional observations. 

3.4. HD 103197: An intermediate mass planet 

We have gathered 58 high quality radial velocity measurement of 
HD 103197 over a time period of 6.12 years. They are shown in 
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Fig. 3. Top panel: Generalized Lomb Scargle periodogram of the 
radial velocity measurements of HD 90156. A FAP of 0.1% is 
indicated by the dashed-dotted line. The peak at about 50 days 
is very clear and has a FAP much less than this value. Middle 
panel: Radial velocities and residuals as a function of Julian 
Date, and (bottom panel) as a function of orbital phase for the 
one-planet Keplerian fit (P = 49.8 days). 



the top panel of Fig.|4] The mean exposure time for this star was 
shorter than for the other three, as it does not belong to the high 
precision sample, and of the order of 5 minutes, which should 
however still be just sufficient to mean out the most important 
stellar oscillation period for this type of star (cf. Pepe & Lovis, 
2008 ). Corresponding to the smaller exposure time, the spectra 
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Table 4. Orbital and physical parameter of the 1 -planet 
Keplerian solution for the companion around HD 103197. The 
eccentricity of the planet is fixed to zero. The epoch of A is JDB 
54846.0994. 



Parameters 




HD 103197 b 


P 


[days] 


47.84 ± 0.03 


A 


[deg] 


151 ±3 


e 




0.0 (fixed) 


CO 


[deg] 


0.0 (fixed) 


K 


[m/s] 


5.9 ±0.3 


V 


[km/s] 


-4.3282 ± 0.0002 


f{m) 


[10~ 9 M Q ] 


0.0010 ±0.0002 


M sin i 


[M s ] 


31.2 ±2.0 


a 


[AU] 


0.249 ± 0.004 


N 




58 


span 


[years] 


6.12 


lt(0 - C) 


[m/s] 


1.4 


Xr 




1.02 



have lower typical S/N at 550 nm, namely between 33 and 107 
with a mean of 59. The error on a single radial velocity mea- 
surement which includes photon noise, calibration and spectro- 
graph drift uncertainty varies between 0.7 and 2.3 m/s, with a 
mean of 1.1 m/s. The observed raw rms of the radial velocities 
is with 3.95 m/s clearly larger, and indicates the presence of a 
companion. As before, we add quadratically additional errors of 
0.8 and 0.5 m/s. For this star, no ThAr simultaneous reference 
mode was usually used, which is however unproblematic given 
the very small nightly drift of the instrument and the fact that we 
mix observations with and without ThAr during one night. 

A periodogram of our radial velocity data shows a very 
strong peak at a period of about 48 days, as can be seen in the 
top panel of Fig. |4j where horizontal lines again indicate FAPs 
of 10 -1 , 10~ 2 and 10~ 3 from bottom to top. Fitting the data with 
one Keplerian orbit leads to a solution with an insignificant ec- 
centricity (0.04 ± 0.05) so that we fix it to zero, leaving us with 
four free parameters. 

The following solution is then found for the companion, as- 
suming a additional radial velocity jitter of 0.8 m/s: An orbital 
period of 47.84 + 0.03 days, a mean longitude A = 151 + 3 deg 
at JDB 54846.0994 and K = 5.9 ± 0.3 m/s, implying a minimum 
mass Msin/ = 31.2 ± 2.0 M @ and a semi -major axis a - 0.25 
AU. The mass corresponds to 1.8 M%. The fit has a reduced;^ 
of 1.02 and cr(0 - C) = 1.40 m/s, which corresponds well with 
the total assumed mean error of the individual measurements. As 
for the two previously discussed stars, using an additional jitter 
of 0.5 m/s instead of 0.8 m/s has only very minor influences on 
the fit, and results in & x 2 r = 1.18. The middle and lower panel 
of Fig. [4] show the radial velocities and the fit as a function of 
time, and phase folded, as well as the residuals around the fit. 
The orbital elements for the Keplerian fit with e — are given in 
Table |4] The errors were again derived from 5000 Monte Carlo 
simulations. 

We note that the orbital period is similar to the rotation pe- 
riod of the star P 10t which we estimate to be 51 +5 days. We 
however checked for periodicity in the bisector shape at such 
a period, and for a significant correlation of the velocity or the 
residual and the bisector shape, finding no such signs. We finally 
looked at the periodogram of the residual around the fit and dis- 
covered no additional significant signal. 
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Fig. 4. Top panel: Generalized Lomb Scargle periodogram of the 
58 radial velocity measurements of HD 103197. The middle and 
bottom panel show the radial velocity and the residuals as a func- 
tion of Julian Date, and as a function of orbital phase together 
with the orbital solution given in table |4] 



4. Summary 

We have reported in this paper the discovery of three extraso- 
lar planet candidates discovered with the HARPS Echelle spec- 
trograph mounted on the 3.6-m ESO telescope located at La 
Silla Observatory. HD 85390b is an intermediate mass planet 
(Msin/ = 42.0 M m ) on an orbit with a semi-major axis of 
1.5 AU. A drift indicates the presence of an additional outer 
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companion, but we cannot currently determine its nature with 
our data. HD 90156 b is a warm Neptune analogue (Msin/ = 
1.05 M-g, a = 0.25 AU). HD 103197 b is another intermediate 
mass planet (Msin i = 31.2 Ms). Its semi-major axis is 0.25 AU. 



5. Discussion: Formation of intermediate mass 
planets 

With intermediate mass planets we here refer to planets with a 
mass larger than Neptune, but less than Saturn. They therefore 
fall in the significant mass gap of a factor 5 between Uranus 
(17.2 M m ) and Saturn (95.2 M e ) which corresponds to an impor- 
tant difference in internal composition, too: Uranus and Neptune 
are ice giants which consist mainly of a heavy elements (iron, 
rocks and ices) and only about 5-15% in mass of hydrogen and 
helium (e.g. Figueira et al. 120091 ). Saturn and Jupiter are in con- 
trast gas giants with only a small fraction of heavy elements. 
Instead, their massive gas envelopes account for most of the 
mass, namely for about 67 to 80% for Saturn, and 87 to 97% 
for Jupiter (Guillot [T999l . 

This difference in mass and composition is in turn under- 
stood in theoretical formation models based on the core accre- 
tion paradigm as the consequence of an important difference in 
the formation history of these two groups of planets: Jupiter and 
Saturn are thought to have undergone a phase of rapid "runaway" 
gas accretion, during which they accreted their large gaseous en- 
velopes on a relatively short timescale, whereas the ice giants 
never underwent this phase (Pollack et al. 1 1996b . 

If the transformation phase from a Neptunian mass to a 
Jovian mass is very short compared to the lifetime of the pro- 
toplanetary disk, i.e. if the gas accretion rate onto the planet is 
very high, then it is very unlikely that the protoplanetary disk 
disappears just during the phase, which would obviously termi- 
nate the gas accretion and leave the planet at some intermediate 
mass. We thus see that in this scenario, the absence of a planet 
with a mass between Uranus and Saturn is not a surprise, but 
rather an expected outcome (Ida & Lin 20Q3). 

One can therefore say that the higher (smaller) the gas ac- 
cretion rate in runaway, the smaller (higher) the expected rela- 
tive frequency of intermediate mass planets compared to giant 
planets. This correlation will manifest itself in the observed de- 
tection rate of intermediate mass extrasolar planets. Indeed point 
detections made over the last years in particular with HARPS to 
a minimum of the observed planetary mass function at about 30- 
40 Mb (e.g. Bouchy et al. 2009). Such a minimum (if further 
confirmed) would be a very important observational finding, as 
it cannot be explained by an observational bias. This is because 
the lower mass Neptunian and Super-Earth planets are more dif- 
ficult to detect than intermediate mass planets. It is interesting to 
note that a relative paucity of intermediate mass planets seems 
to exist also in transit searches (Hartman et al. 2009). 

Several physical effects influence the runaway gas accretion 
rate (Lubow et al.[T999l Ida & Lin|2003 Mordasini et al. |2009al 
Lissauer et al. 2009): Thermal pressure in the planetary enve- 
lope at lower planetary masses, the global evolution (dissipation) 
of the protoplanetary disk which reduces the amount of matter 
available for the planet, the rate of mass transport due to viscos- 
ity within the disk towards the planet as well as local phenomena 
like the formation of a gap. 

While several planetary population synthesis simulations 
built on the core accretion paradigm share the common feature 
that they predict some depletion of intermediate mass planets 
relative to other planet types, they attribute in their underly- 




100 1000 
Msini [M e ] 

Fig. 5. Theoretical mass distribution from Neptunian to Jovian 
mass planets obtained from population synthesis calculations. 
The solid line shows a population where the planetary gas accre- 
tion rate was only limited by the disk accretion rate if the planet 
has a mass larger than the local gas isolation mass. For the dot- 
ted line, the limit was used in any case. Both distributions were 
normalized to unity at the first bin at about 20 Earth masses. 



ing formation models (cf. Alibert et al. 2005) distinct impor- 
tance to the various effects mentioned above, so that there are 
clear differences in the degree of depletion: from an almost 
complete absence of such planets between 0.1-1 AU (Miguel 
& Brunini 2008 2009 ), over a significant depletion (Ida & Lin 
1200411200 8a. 2008b ) to a rather moderate one (about a factor 2-3 
relative to Jovian planets) in Mordasini et al. (2009a 2009b). 

The reason for the difference between the latter two models 
was discussed in details in Mordasini et al. (2009a). It is related 
to different assumptions about the amount of gas present in the 
vicinity of the planet as discussed next. Here we illustrate this 
in Fig.[5]that shows the theoretically obtained mass distribution 
from Neptunian to Jovian planets using two different assump- 
tions for the gas accretion rate in runaway. Other settings are 
similar as in Mordasini et al. (2009a). Only planets with a pe- 
riod less than 5 years are shown which are detectable for a RV 
instrument of 1 m/s precision, similar as HARPS. 

In one simulation (solid line), the gas accretion rate of the 
planet is limited by the accretion rate in the disk only if the 
mass of the planet is larger than the local gas isolation mass, 
calculated with the undisturbed gas surface density. This is sim- 
ilar to the criterion used by Ida & Lin (2004). It means that 
for masses between roughly 30 up to ~ 100 M @ i.e. before the 
gas isolation mass is reached, high gas accretion rates of up to 
roughly 10~ 2 M e /yr can occur, so that the transformation from 
a Neptunian to a Jovian planet only takes of order 10 4 years, 
much less than typical disk lifetimes. The underlying assump- 
tion is here that gas already inside the planet's Hill sphere can 
be accreted independently of the inflow from further away in the 
disk. This inflow is in turn limited by the disk viscosity. 

In the other simulation (dotted line), the planetary gas accre- 
tion rate is limited by the accretion rate in the disk as soon as 
runaway starts. As cores reach a mass large enough to trigger 
gas runaway accretion at a moment typically not much before 
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the disk goes away (it is recalled that around about 90% of FGK 
stars, there are no giant planets), the accretion rate in the disk has 
usually already fallen to quite low values at this moment, of or- 
der a few 10~ 4 to 1CT 3 M e /yr. This means that the transformation 
from a Neptunian to a Jovian planet takes now several 10 5 years, 
which is no more much shorter than the remaining disk life time. 
The underlying assumption is here that due to gap formation, the 
mass directly available to the planet is in fact small, in particular 
smaller than the gas isolation mass as calculated above, as (be- 
ginning) gap formation reduces the gas surface density around 
the planet. This corresponds to the setting in Mordasini et al. 
(20091Q. 

The plot shows that the two different settings have an impor- 
tant effect on the frequency of intermediate mass planets. They 
are in the first model clearly less frequent that in the second one, 
with a factor of two of difference for M * 40M ffi . Thus, the 
frequency of intermediate mass planets can be used to directly 
improve our theoretical understanding of the accretion process. 

Bodies with masses of up to several ten Earth masses can 
in principle also form after the dissipation of the gas disk if 
enough solids are available in situ, which is in particular the case 
at larger distances beyond the ice line. Such bodies would then 
be essentially gas free. At the rather small orbital distances of 
HD 85390b and HD 103197b (1.5 and 0.25 AU, respectively), 
for realistic solid disk masses, such a formation scenario seems 
however unlikely: The two planets are larger than what can be 
formed by this process at their current distances (Ida & Lin 
2004), considering the host star metallicities. 

The two intermediate mass planets were therefore very prob- 
ably formed while the gas disk was still present. In this case, 
the mass of a planet core cannot grow to arbitrarily large val- 
ues without runaway gas accretion setting in. The real masses of 
HD 85390 b and HD 103197 b which are on the statistical mean 
a factor 4/7T larger than the projected ones are probably larger 
than the mass at which this process starts (although not by a 
large factor, especially for HD 103197 b, and the specific value 
depends on e.g. the unknown core accretion rate, see Papaloizou 
& Terquem 119991 1. The two planets are therefore probably ex- 
amples where gas runaway started, but only shortly before the 
gas disk disappeared, so that only low quantities of gas were 
still available to accrete, and that the gas accretion rate was low. 
Indeed consist synthetic planets (Mordasini et al. 2009b| that 
are situated at a similar position in the mass-distance plane as 
HD 85390 b (M sin i = 42 M e ) of typically about 40% hydrogen 
and helium in mass. The scatter around this value is quite large, 
reflecting different disk properties, and ranges from still clearly 
solid dominated planets (~ 10% gas) to small gas giants (~ 70% 
gas). For the less massive HD 103197 b (Msini = 31.2 M e ) the 
typical value is about 30% of gas, with a scatter around this 
value of about 15%. We conclude that these two exoplanets have 
thus probably not only a mass, but also a composition between 
Neptunian and Jovian planets. 

From this discussion we see that it is important to observa- 
tionally infer relative frequencies of planetary types between 20 
to 100 M m to better understand the runaway phase. This will fi- 
nally enable us to construct better formation models. In the ideal 
case, the planets should be transiting their host star, but be lo- 
cated still so far away from it that significant evaporation can be 
ruled out, so that the primordial composition can be studied. 
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